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Abstract 

We present the results of experimental and numerical study of the distribution of the reflection 
coefficient P(R) and the distributions of the imaginary P(v) and the real P{u) parts of the Wigner's 
reaction K matrix for irregular fully connected hexagon networks (graphs) in the presence of strong 
absorption. In the experiment we used microwave networks, which were built of coaxial cables 
and attenuators connected by joints. In the numerical calculations experimental networks were 
described by quantum fully connected hexagon graphs. The presence of absorption introduced 
by attenuators was modelled by optical potentials. The distribution of the reflection coefficient 
P(R) and the distributions of the reaction K matrix were obtained from the measurements and 
numerical calculations of the scattering matrix S of the networks and graphs, respectively. We 
show that the experimental and numerical results are in good agreement with the exact analytic 
ones obtained within the framework of random matrix theory (RMT). 
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Quantum graphs of connected one-dimensional wires were introduced more than sixty 

n n 

years ago by Pauling pj. Next the same idea was used by Kuhn [2J to describe organic 
molecules by free electron models. Quantum graphs can be considered as idealizations of 
physical networks in the limit where the lengths of the wires are much bigger than their 
widths, i.e. assuming that the propagating waves remain in a single transversal mode. 
Among the systems modelled by quantum graphs one can find e.g., electromagnetic optical 
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soscopi 



quantum wires 



and excitation of fr actons 



in fractal structures 
paradigms of quantum chaos 



Recently it 
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ias been shown that quantum graphs are excellent 
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2l|. More complicated 



and thus more realistic systems - microwave networks with moderate absorption strength 
7 = 2nr I A < 7.1, where T is the absorption width and A is the mean level spacing, have 



been experimentally investigated in 22j, |23j. Other interesting open objects - quantum 
graphs with leads - have been analyzed in details in Q, Q. However, the properties of 
networks and graphs with strong absorption have not been studied experimentally neither 
numerically so far. Therefore, in this paper we study experimentally and numerically the 
distribution of the reflection coefficient P{R) and the distributions of the Wigner's reaction 



matrix 



24j (in the literature often called K matrix [25|]) for networks (graphs) with time 



reversal symmetry (f3 = 1) in the presence of strong absorption. 

In the case of a single channel antenna experiment the K matrix is related to the scattering 
matrix S by the following relation 

S=l^. (1) 



1 + iK 



Eq. (1) holds for the systems with absorption but without direct processes 25] . It is 
important to mention that the function Z = iK has a direct physical meaning of the electric 
impedance that has been recently measured in the microwave cavity experiment |26(]. In the 
one channel case the S matrix can be parameterized as 



S 



Re 



w 



(2) 



where R is the reflection coefficient and 9 the phase. 

Properties of the statistical distributions of the S matrix with direct processes and im 



per 
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5ect coupling have been studied theoretically in several important papers 27 
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30, 



321 ] . Recently the distribution of the S matrix has been also measured experimentally 



for chaotic microwave cavities with absorption 



331 ] . The distribution P(R) of the reflec- 



tion coefficient R and the distributions of the imaginary P(v) and the real P(u) parts of 
the Winner's reaction K matrix are theoretically known for any dimensionless absorption 
3^ |. In the case of time reversal systems (symmetry index (3 = 1) P(R) 



strength 7 



25 



has been studied experimentally by Mendez- Sanchez et al 
and P(u) have been studied for chaotic microwave cavities in 
networks for moderate absorption strength 7 < 7.1 in 



35 L The distributions P( 
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26l . l36| and for microwave 



231 ] . For systems without time 



reversal symmetry (/3 = 2) and a single perfectly coupled channel P(R) was calculated by 



Beenakker and Brouwer 



371 ] while the exact formulas for the distributions P(v) and P(u) 



were given by Fyodorov and Savin [251 ]. 

In the experiment quantum graphs can be simulated by microwave networks. The analogy 
between quantum graphs and microwave networks is based upon the equivalency of the 
Schrodinger equation describing the quantum system and the telegraph equation describing 



the microwave circuit 



2l|. 



A general microwave network consists of N vertices connected by bonds e.g., coaxial 
cables. A coaxial cable consists of an inner conductor of radius r\ surrounded by a concentric 
conductor of inner radius Ti. The space between the inner and the outer conductors is filled 
with a homogeneous material having a dielectric constant e. For a frequency v below the 
onset of the next TEn mode only the fundamental TEM mode can propagate inside a coaxial 
cable. (This mode is in the literature often called a Lecher wave.) The cut-off frequency 



of the TE11 mode is v c 



n(ri+r2)\/e 



32.9 GHz 



where r% = 0.05 cm is the inner wire 



radius of the coaxial cable (SMA-RG402), while T2 = 0.15 cm is the inner radius of the 
surrounding conductor, and e ~ 2.08 is the Teflon dielectric constant {39} Q]. 

From the experimental point of view absorption of the networks can be changed by the 
change of the bonds' (cables') lengths [21I or more effectively by the application of microwave 



attenuators 
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23] . In the numerical calculations weak absorption inside the cables can be 



211 ] . We will show that strong absorption 



described with the help of complex wave vector 
inside an attenuator can be described by a simple optical potential. The corresponding 



mathematical theory has been developed in 411 ] . 



The distribution P(R) of the reflection coefficient R and the distributions of the imaginary 
and real parts of the Wigner's reaction matrix K for microwave networks with absorption 



were found using the impedance approach 
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26 



36] . In this approach the real and imagi 
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FIG. 1: (a) The scheme of the experimental set-up for measurements of the scattering matrix S n 
of the microwave fully connected networks with absorption. Absorption in the networks was varied 
by the change of the attenuators, (b) The scheme of the setup used to measure the radiation 
scattering matrix S r . Instead of a network five 50 f2 loads were connected to the 6-joint. 



nary parts of the normalized impedance Z 

Re Z n + i(lm Z n - Im Z r ) 

of a chaotic microwave system are measured, with Z n t r ) = Z (l + S n ^)/ (1 — S n ( r )) being the 
network (radiation) impedance expressed by the network (radiation) scattering matrix S n t r ) 
and Zq is the characteristic impedance of the transmission line. The radiation impedance 
Z r is the impedance seen at the output of the coupling structure for the same coupling 
geometry, but with the vertices of the network removed to infinity. The Wigner's reaction 
matrix K can be expressed by the normalized impedance as K = —iZ. The scattering 
matrix S of a network for the perfect coupling case (no direct processes present) required 
for the calculation of the reflection coefficient R (see Eq. (2)) can be finally extracted from 
the formula S = (1 - Z)/(l + Z). 

Figure 1(a) shows the experimental setup for measuring the single-channel scattering ma- 
trix S n of fully connected hexagon microwave networks necessary for finding of the impedance 
Z n . We used Hewlett-Packard 8720A microwave vector network analyzer to measure the 
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scattering matrix S n of the networks in the frequency window: 7.5-11.5 GHz. The networks 
were connected to the vector network analyzer through a lead - a HP 85131-60012 flexible 
microwave cable - connected to a 6-joint vertex. The other five vertices of the networks were 
connected by 5-joints. Each bond of the network presented in Fig. 1(a) contains a microwave 
attenuator. 

The radiation impedance Z r was found experimentally by measuring the scattering matrix 
S r of the 6-joint connector with five joints terminated by 50 Q loads (see Figure 1(b)). 

The experimentally measured fully connected hexagon networks were described in nu- 
merical calculations by quantum fully connected hexagon graphs with one lead attached 
to the 6-joint vertex. In the calculations attenuators (absorbers) were modelled by optical 



potentials 41] . To be explicit we suppose that the fully connected hexagon graph T with 
one coupled antenna is described in the Hilbert space L 2 (T) := 0Q >n ) L 2 (0, £j n ) L 2 (0, oo), 
where £j n stays for the lengths of the bond connecting the vertices j and n and the halfline 
(0, oo) describes the attached antenna. 

We define the Schrodinger operator H by 

Hif; jn := -ip'! n + U jn ip jn , (4) 

with ipj n G L 2 (0,£j n ) for the bonds and 



H*p 0n := -C (5) 

with ipQ n G L 2 (0, oo) describing the wave function of the antenna connected to the vertex 
n (note that the "infinite" vertex of the antenna has index 0) . 

At the vertices the wave functions are linked together with the boundary values 

iijnti) ■= jim_ ^jn(x) , ip' jn (j) : lim + ip' jn (x) , (6) 

satisfying boundary conditions 4>j n (j) = 4>jm(j) ='■ 4>j for all n,m describing connected 
vertices, and 

E ^nU) = 0. (7) 
7161/ (5) 

The optical potentials Uj n are purely imaginary and describe the absorber inserted be- 
tween the vertices (j, n). 
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Since the graph T is infinite (due to the attached antenna) we can look for solutions of 
the equation 



Hit = k 2 tb, (8) 

referring to the continuous spectrum, where k is the wave vector. For microwaves propa- 
gating inside a lossless bond with a dielectric constant e the wave vector k = 2irei>/c, where 
v and c denote the frequency of a microwave field and the speed of light in the vacuum, 



respectively. To solve this equation we used the graph duality principle [4l|. According 
to this principle we need to solve the equation — /" + Uj n f = k 2 f on [0,€j n ] satisfying the 
normalized Dirichlet boundary conditions 

u jn (£ jn ) = l-(u jn )'(l jn ) = , v jn (0)l-(v jn )'(0) = . (9) 
The Wronskian of this solution is naturally equal to Wj n — Vj n (£j n ) = Uj n (0). Then 



according to [41[ the corresponding boundary values (J6J) satisfy the equation 



^" r l* = ' (10) 



Conversely, any solution ipj of the system (|T0|) determines a solution of (jHJ) by 

IPjn(x) = ~tt7~~ Ujn(x) - V jn (x) if U = 1, .., 6 , (11) 

lb- 

ib jn (x) = - — — v jn {x) if n = . (12) 

Wj n 

As already mentioned the microwave attenuators are modelled by optical potentials local- 
ized inside the inserted component. It is well known that any smooth and localized potential 
can be easily approximated by a sequence of delta potentials inside the support of the poten- 
tial - see ||42j for details. We will use this fact and express the optical potential as a sum of N 
delta-potentials with imaginary coupling constants: U(x) = ibJ2r=i <K X ~~ ( r — l)h/{N — 1)). 
The delta-potentials were equally spaced inside the length If, of the absorbing element (at- 
tenuator). By changing the number N and the strength b of delta-potentials we were able 
to vary absorbing properties as well as reflective properties of attenuators. We used N = 10 
delta-potentials with b = 0.028 m -1 for simulation of the 1 dB attenuators and N = 12 
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FIG. 2: In the panels (a) and (b) the modulus \S n \ and the phase 6 of the scattering matrix S 
measured for the network with 7 = 19.9 are plotted in the frequency range 7.5 - 9 GHz. In (c) 
and (d) \S n \ and 9 of the scattering matrix S n are plotted for the network with 7 = 47.9 in the 
same frequency range. The measurements have been done for the two networks which in each bond 
contained: 1 dB attenuator ((a) and (b)) and 2 dB attenuator ((c) and (d)), respectively. The 
total "optical" length of the microwave networks including joints and attenuators were 574 cm and 
554 cm, respectively. 



delta-potentials with b = 0.045 vrC 1 for the 2 dB attenuators. In both cases the length of 
the attenuator was = 2.65 cm. Furthermore, in the numerical calculations of the scatter- 
ing matrices S n of the graphs the weak absorption inside the microwave cables was taken 
into account by replacing the real wave vector k by the complex vector k + iay/k Q|, where 
the absorption coefficient was assumed to be a = 0.009 m -1 / 2 21], 

In order to find the distribution P(R) of the reflection coefficient R and the distributions of 
the imaginary and real parts of the K matrix we measured the scattering matrix S n of 88 and 
74 network configurations containing in each bond a single 1 dB and 2 dB microwave SMA 
attenuator, respectively. The total optical lengths of the microwave networks containing 
1 dB attenuators, including joints and attenuators, varied from 574 cm to 656 cm. For 
the networks with 2 dB attenuators the optical lengths varied from 554 cm to 636 cm. To 
avoid degeneracy of eigenvalues of the networks the lengths of the bonds were chosen as 
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incommensurable . 

In Figure 2 the modulus \S n \ and the phase 6 of the scattering matrix S n of the microwave 
networks with 7 = 19.9 and 47.9, respectively, are presented in the frequency range 7.5 - 9 
GHz. The measurements were done for two networks containing 1 dB and 2 dB attenuators, 
respectively. Their total "optical" lengths including joints and attenuators were 574 cm and 
554 cm, respectively. 

For systems with time reversal symmetry ((3 = 1), the explicit analytic expression for the 
distribution P(R) of the reflection coefficient R is given by [34] 



The probability distribution Pq(x) is given by the expression 



(13) 



Po(x) 



dW{x) 
dx 



(14) 



where the integrated probability distribution W(x) is expressed by the formula 34 ] 



W{x) 



X + 1 r 



47T 



fx(w)g 2 (w) + f 2 (w)g 1 (w) + h 1 {w)j 2 {w) + h 2 (w)j 1 (w) 



w=(x-l)/2 



The functions fi,g%, h\,ji are defined as follows 



(15) 



dt 



t\t-w \e-~< t/2 

(1 + t)3/ 2 



1 - e" 7 + - 



(16) 



hAw) 



dt 



dt 



-7*/2 



t \ t-W |(l + t) 3 / 2 

^\t-w\e~^ 2 



y/t(l + 1) 



7+(l-e-^)( 7 t-2) 



3i{w) 



dt- 



-7*/2 



(17) 



(18) 



(19) 



\Jt\t- w \(i + ty/ 2 ' 

Their counterparts with the index 2 are given by the same expressions but the integration 
is performed in the interval t G [0, w] instead of [w, 00). 



S 




FIG. 3: Experimental distribution P(R) of the reflection coefficient R for the microwave fully 
connected hexagon networks at 7 = 19.3 (open squares) and 7 = 47.7 (full squares). The cor- 
responding theoretical distribution P(R) evaluated from the Eq. (|13p is marked by the solid line 
(7 = 19.3) and dashed line (7 = 47.7), respectively. 



DC 




FIG. 4: Numerical distribution P(R) of the reflection coefficient R for fully connected hexagon 
quantum graphs at 7 = 19.3 (open circles) and 7 = 47.7 (full circles). The corresponding theoretical 
distribution P(R) evaluated from the Eq. (I13D is marked by the solid line (7 = 19.3) and dashed 
line (7 = 47.7), respectively. 



9 



The distributions of the imaginary and the real parts P(y) and P(u) of the K matrix 



25( can be also expressed by the probability distribution Pq(x): 



V2 



7TV- 

and 

1 



1 / 1 



^s^f^F^V;)]- < 21) 

where — v = ImK < and u = KeK are, respectively, the imaginary and real parts of the 
K matrix. 

Figure 3 shows the experimental distributions P(R) (squares) of the reflection coefficient 
R for two mean values of the parameter 7, viz., 19.3 and 47.7. The distribution for 7 = 19.3 
is obtained by averaging over 88 realizations of the microwave networks containing 1 dB 
attenuators. The distribution for 7 — 47.7 is obtained by averaging over 74 realizations 
of the microwave networks containing 2 dB attenuators. The experimental values of the 7 
parameter were estimated for each realization of the network by adjusting the theoretical 
mean reflection coefficient (R)th to the experimental one (R) = (SS^), where 

(R) th = f 1 dRRP(R). (22) 
Jo 

Figure 3 also presents the corresponding distributions P(R) (solid and dashed lines, 
respectively) evaluated from Eq. ({TBI . A good overall agreement of the experimental distri- 
butions P(R) with their theoretical counterparts is seen. 

Figure 4 shows the numerically evaluated distributions P(R) (circles) of the reflection 
coefficient R for the graphs at 7 = 19.3 and 47.7 compared to the theoretical ones evaluated 
from the formula Eq. (fl3l) . The numerical distributions are the result of averaging over 
162 and 214 realizations of the graphs with optical potentials simulating 1 dB and 2 dB 
attenuators, respectively. The numerical values of 7 parameter were also estimated by 
adjusting the theoretical mean reflection coefficient to the numerical one. The agreement 
between the numerical results for 7 = 47.7 and the theoretical ones (dashed line) is good. 
However, for 7 = 19.3 for R < 0.15 some discrepancies between the numerical results and 
the theoretical ones (solid line) are visible. 

In Figure 5 the experimental distribution P(v) of the imaginary part of the K matrix 
is shown for the two mean values of the parameter 7 = 19.3 and 47.7, respectively. The 
distribution is the result of averaging over 88 and 74 realizations of the networks with the 
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FIG. 5: Experimental distribution P(v) of the imaginary part of the K matrix for the two values of 
the mean absorption parameter: 7 = 19.3 (open squares) and 7 = 47.7 (full squares), respectively. 
The corresponding theoretical distribution P(v) evaluated from the Eq. (|20p is marked by the solid 
line (7 = 19.3) and dashed line (7 = 47.7), respectively. 




FIG. 6: Numerical distribution P(y) of the imaginary part of the K matrix for the two values of 
the mean absorption parameter: 7 = 19.3 (open circles) and 7 = 47.7 (full circles), respectively. 
The corresponding theoretical distribution P(v) evaluated from the Eq. (120ft is marked by the solid 
line (7 = 19.3) and dashed line (7 = 47.7), respectively. 
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attenuators 1 dB and 2 dB, respectively. The experimental results in Figure 5 are in general 
in good agreement with the theoretical ones. However, both experimental distributions are 
slightly higher than the theoretical ones in the vicinity of their maxima. 

Results of the numerical calculations of the distributions P(v) are shown in Figure 6 
for two mean values of the parameter 7 = 19.3 and 47.7, respectively. They are compared 
to the theoretical ones evaluated from the formula Eq. (|20|) . Figure 6 shows also a good 
agreement between the numerical and theoretical results, which confirms usefulness of the 
optical potential approach in describing the microwave attenuators. 

Measurements of the distribution P(u) of the real part of the Wigner's reaction matrix 
give an additional test of the consistency of the 7 evaluation. In Figure 7 we show this 
distribution obtained for two values of 7 = 19-3 and 47.7, respectively, compared to the 
theoretical ones evaluated from the formula Eq. (|2ip . Also here we observe good overall 
agreement between the experimental and theoretical results. However, Figure 7 shows that 
for the networks with 2 dB attenuators the theoretical distribution is in the middle (—0.1 < 
u < 0.1) slightly higher than the experimental one. According to the definition of the 
K matrix (see Eq. (1)) such a behavior of the experimental distribution P(u) suggests 
deficiency of small values of |ImS|, whose origin is not known. Moreover, the experimental 
distribution P(u) obtained for the networks assembled with 1 dB attenuators is slightly 
asymmetric for \u\ > 0.5. 

In Figure 8 the comparison of the numerical distribution P(u) obtained for two values of 
7 = 19.3 and 47.7, respectively, to the theoretical one evaluated from the formula Eq. ( |2~T1) 
is presented. In this case we see a good overall agreement between the numerical and 
theoretical results. 

In spite of the above mentioned discrepancies which appeared mainly in the case of the 
experimental distribution P(u) the overall good agreement between the experimental and 
theoretical results justifies a posteriori the chosen procedure of calculating the experimental 
7. The same is true also for the numerical simulations. 

In summary, we measured and calculated numerically the distribution of the reflection 
coefficient P(R) and the distributions of the imaginary P(v) and the real P{u) parts of 
the Wigner's reaction K matrix for irregular fully connected hexagon networks and graphs 
in the presence of strong absorption. In the case of the microwave networks consisting of 
SMA cables and attenuators the application of attenuators allowed for effective change of 



12 



1 .6 i 1 1 1 1 ' 1 ' 1 1 r 




U 

FIG. 7: Experimental distribution P(u) of the real part of the K matrix for the two values of the 
mean absorption parameter: 7 = 19.3 (open squares) and 7 = 47.7 (full squares), respectively. The 
experiment is compared to the theoretical distribution P{u) evaluated from the Eq. (j 2 1|) : solid line 
(7 = 19.3) and dashed line (7 = 47.7). 




FIG. 8: Numerical distribution P(u) of the real part of the K matrix for the two values of the 
mean absorption parameter: 7 = 19.3 (open circles) and 7 = 47.7 (full circles), respectively. The 
numerical results are compared to the theoretical distributions P{u) evaluated from the Eq. (I2ip : 
solid line (7 = 19.3) and dashed line (7 = 47.7). 
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absorption in the graphs. In the numerical calculations absorption in an attenuator was 
modelled by an optical potential. We showed that in the case of the time reversal symmetry 
(f3 = 1) the experimental and numerical results for P(R), P(v) and P(u) are in good overall 
agreement with the theoretical predictions. The agreement of the numerical and theoretical 
results strongly confirms the usefulness of the optical potential approach in the description 
of the microwave attenuators. 
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